A Candidate Brightest Proto-Cluster Galaxy at z = 3.03 
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ABSTRACT 

We report the discovery of a very bright (hir = 22.2) Lyman break galaxy at 
z = 3.03 that appears to be a massive system in a late stage of merging. Deep 
imaging reveals multiple peaks in the brightness profile with angular separations 
of ~0."8 (~25 h^ 1 kpc comoving). In addition, high signal-to-noise ratio rest- 
frame UV spectroscopy shows evidence for ~ 5 components based on stellar 
photospheric and ISM absorption lines with a velocity dispersion of a ~ 460 km 
s _1 for the three strongest components. Both the dynamics and high luminosity, 
as well as our analysis of a ACDM numerical simulation, suggest a very massive 
system with halo mass M ~ 10 13 M Q . The simulation finds that all halos at z = 3 
of this mass contain sub-halos in agreement with the properties of these observed 
components and that such systems typically evolve into M ~ 1O 14 M halos in 
groups and clusters by z = 0. This discovery provides a rare opportunity to 
study the properties and individual components of z ~ 3 systems that are likely 
to be the progenitors to brightest cluster galaxies. 

Subject headings: galaxies: high redshift — galaxies: mergers — galaxies: evolu- 
tion — galaxy clusters: high redshift clusters 
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Introduction 



Massive galaxies and galaxy clusters at high redshift provide strong constraints on cos- 
mological models and galaxy formation scenarios. The Lyman break galaxies (LBGs) are 
a widespread population of high redshift galaxies selected by their rest-frame FUV colors 



dsteidel et al.lll99djljjp l LBGs at z ~ 3 tend to reside in massive halos (M DM ) ~ lO n ' 5 M 



jAdelbereer et alJl2005| ) and are therefore believed to evolve into present-day massive ellip 
tical galaxies. Although ~ 3000 color-selected z > 2 LB Gs have spectroscopic confirmations 



;e.g. iSteidel et al J 12003 . 12004 : ICooke et al.ll2005l . 12008 ). only two systems, selected by their 
radio emission, have been spectrosc opically identified as potential brightest cluster galaxies 
rtzirm et alihood : iMiley et al-lbood l. 



We report the discovery of a bright, m# = 22.2, rest-frame FUV color-selected LBG 
(hereafter LBG-2377) at z = 3.03. Ground-based imaging and spectroscopy finds no evi- 
dence for AGN activity or gravitational lensing making this the most luminous z > 3 LBG 
discovered to date. The deep images and high signal-to-noise ratio (SNR) spectroscopy 
show evidence for multiple massive components indicative of a system in the late stage of 
merging. We show that the component velocity dispersion and integrated FUV luminosity 
suggest that LBG-2377 resides in a very massive halo (~ 1O 13 M ) and is likely to evolve into 
a present-day massive brightest cluster galaxy (§ 13. 2p . Therefore, LBG-2377 provides a rare 
opportunity to study the dynamics and properties of a massive LBG halo in the process of 
formation. 

All z > 3 LBGs with R < 23 discovered to date are magnified by gravitational lensing. 
Because all of these show Ly-a in absorption, they are representative o f the spectral proper- 
ties observed in ~25-50% of the LBG population (jShapley et al.ll2003l ). LBG-2377 displays 
strong Ly-a in emission in two of its components and, because of its brightness, is the only 
system of its kind. The luminosity of LBG-2377 enables high SNR study of its individual 
components and thereby a means to constrain the observed range of LBG spectral properties. 
We briefly describe the imaging and spectroscopic observations of LBG-2377 in §HJ highlight 
the initi al analysis of the dat a in §31 and present a discussion in §U All magnitudes are in 
the AB ( IFukugita et al.lll996l ) magnitude system unless otherwise noted and we assume an 
fijif = 0.3, = 0.7 cosmology and use comoving units throughout. 
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Observations 



2.1. Imaging 



The deep images of the QSO PC1643+ 4631A field were obtain ed on 2001 April 18 using 
the Low- Resolution Imaging Spectrometer (ILRIS; Oke et al.lll995l ) mounted on the Keck I 
telescope as part of a multi- field survey for LBGs associated with damped Ly-ct systems 
at z ~ 3 (ICooke et al.ll2005l ). Specifically, the field containing LBG-2377 was imaged in 
the u filters with seeing of 0.7 — 0.9 arcsec and a depth of u' = 28.8 and BVRI 

> 28.1 (la). We discovered LBG-2377 as an R = 22.2 z ~ 3 LBG candidate at R.A. 
16:44:48.3, Dec.+46:27:08.2 (J2000.0) towar d one edge of the image d field. The isophotal 
magnitudes as determined using S Extractor (IBertin fc Arnoutslll996l ) source extraction code 
are: v! = 26.5 ± 0.46, B = 24.3 ± 0.10, V = 22.6 ± 0.03, R = 22.2 ± 0.06, and I = 22.3 ± 0.07. 



Gardner et al.ll200ll ). As illus- 



LBGs at z ~ 3 have a half-light radius of < 0."3 (e.g., 
trated in the contour plot shown in Figure [fl LBG-2377 spans ~ 2" x 3" and has a morphology 
that can be interpreted as a single extended system exhibiting multiple vigorous star forming 
regions or multiple separate near-point sources. The ratio of the peaks of the three strongest 
identified components is approximately 10:5:1 with separations of ~0."6 — 1".0 (~ 20 — 30ft -1 
kpc, comoving). Inspection of the deep images shows that there are no significant lensing- 
source candidates within ~ 10". A potential candidate 20" from LBG-2377 (R.A. 16:44:46.5, 
Dec +46:27:04.1) does have an effective redshift for lensing (z p hat = 0.41 ± 0.2), but would 
need a halo m ass of M > 7 x 1O 14 M to provide a significant (> 2 — 3x) magnification boost 
( lMinorll2007l ). Moreover, such a geometry would result in a large distortion of the LBG-2377 
image. We see no evidence for this in the deep ground-based images. 



2.2. Spectroscopy 



LBG-2377 was confirmed to be a z = 3.03 galaxy from l ow-resolution discovery multi- 
object spectroscopy on 2004 February 18 (ICooke et al.ll2005l ). The data were taken using 
the 300 line mm -1 grism on LRIS with a spectral resolution of ~ 10A FWHM. We acquired 
follow-up higher resolution, higher SNR longslit spectroscopic observations with LRIS on 21 
May 2007 using the 600 line mm -1 grism blazed at 4000A with the blue arm and the 1200 



1 R S is the H. Spinrad high-throughout night-sky R filter publicly available at the W. M. Keck Observatory. 
The tr ansformation from Spinrad night-sky R s into Johnson VR is (R s - R) = -0.004 - 0.072(V - R) + 0.073(V 



R) 2 ( Stern et al.lll999h . All Rs magnitudes presented here have been transformed into the Johnson R 



bandpass. 
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line mm 1 blazed at 7500 A with the red arm. We used five 1800s exposures in our final 
combined spectrum that resulted in a SNR of 10- 15 and a resolution of ~ 2-3A(Figure[2]). 
Details of the observations will appear in a forthcoming paper. 

The spectrum of LBG-2377 is dominated by O and B star continua and displays strong 
Ly-a emission and little reddening. These features are cons istent with the general properties 



of the z ~ 3 LBG population that display Ly-a in emission (iShapley et al.ll2003l : ICooke et al. 



20051 ). Although the Ly-a feature appears as a single emission line in the low- resolution 
discovery spectrum, the higher-resolution high SNR spectroscopy reveals two strong peaks 
with two possible weak peaks more evident in the individual spectra. The rest-frame FUV 
shows a complex series of interstellar atomic lines exhibiting gas absorption over velocities 
> 2000 km s _1 and evidence for ~ 5 components. Over 20 FUV ISM absorption lines were 
used to identify and study individual components. 



3. Analysis 
3.1. Components 



We identified >5 photospheric lines (to varying confidence levels) of each component to 
estimate their systemic redshifts. These lines include^: Sim A1294.54, 1296.73, 1417.24, Cm 
A1175.26, 1296.33, 1427.85 Cn A1323.93, Nm A1324.32, Oiv A1343.35, and Sv A1501.76. 
The strongest components have redshifts z = 3.0385, 3.0289, and 3.0244. We find a 1-D 
velocity dispersion of a = 460 km s _1 for these three components and a = 705 km s _1 for 
the five components. The details of the identified components are listed in Table [D 

We note that the Ly-a emission and ISM absorption lines of the LBG population show 
observed vel ocity offsets of ~ +650 and ~ —200 km sec" 1 from their sys temic velocities, 



respectively (IPettini et al.l 120021 ; lAdelberger et al.ll2003l ; ICooke et al.ll2005l ). The Ly-a and 
ISM offsets in the identified components are consistent with that of the LBG population. 
The accepted interpretation is that these offsets are caused by the presence of an expanding 
galactic-scale shell of gas and dust driven by supernovae and stellar winds. Ly-a photons are 
largely absorbed in the approaching shell but are resonantly scattered off the receding portion 
and survive. In this scenario, the two prominent Ly-a peaks observed in the spectrum of 
LBG-2377 would trace two distinct expanding galactic-scale shells with a systemic velocity 
difference of ~ 550 km sec -1 . Our fit to the components in Table [1] is shown in Figure [3j 



2 For consistency with future analyses and because our observ ations are of similar quality, we use the 
stellar photospheric values listed in Table 1 of IPettini etall |200d) . 
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Although the spectra of the components are superposed, the ISM absorption-line strengths 
and behavior of each component appear internally consistent with being separate systems. 



3.2. Theoretical expectations 



The high luminosity and complex structure of LBG-2377 suggest that this system may 
be a brightest cluster galaxy in formation. Here we investigate this possibility in the context 
of ACDM. We estimate the dark matter halo mass of LBG-2377 using both its measured 
kinematics and luminosity. Encouragingly, both estimates suggest a group-scale mass ~ 
1O 13 M . If we assume an NFW halo, our measured a ~ 460 km s _1 1-D velocity dispersion 
corresponds to a ha lo maximum circular velocity of V max ~ 1.5a, or V mK * ~ 700 km s" 1 
fclvpin et alJll999h . At z ~ 3 this implies M vir ~ 10 13 M (e.g. iBullock et aJboOlh . 



Alternatively, if we assume a simple monotonic relationship between galaxy luminosity and 
dark matter halo L ,, (e.g. iBerrier et al.l l2006r) and use the z ~ 3 LBG luminosity function 



(ISteidel et al.l Il999l ; ISawicki fc Thompson! 120061 ). the number density of objects similar to 
LBG-2377 (hir = 22.2) is ~ 10~ 7 Mpc~ 3 /t 3 . This n umber density corresponds to halos 
with Vmax ~ 700 km s -1 at z = 3 (IBerrier et al.l 120061 ) . The agreement between these two 
estimates suggests that LBG-2377 does indeed sit within a massive V max ~ 700 km s^ 1 halo. 

Halos that are this massive at z = 3 almost certainly contain substructure. We illustrate 
this explicitl y by examining the high-resolution 80 h~ 3 Mpc~ 3 ACDM N-body simulation 
described in IStewart et al.l (120071 ). There are three halos in this box with V majX > 600 km 
s _1 at z = 3 (specifically, V max = 620,645, and 740 km sec -1 ). These objects have between 
five and eight lO n /i _1 M < M < 10 12,5 /i _1 M Q sub-halos within their (physical) virial radii 
i? vir ~ 90k- 1 kpc (M/lO 13 /i~ 1 M ) 1 / 3 . Note that this substructure mass is ty pical of that 
expected for LBG hosts at z ~ 3 ((M LBG ) ~1O 1L5 M kdelberger et al.ll2005h . If we look 
at simulated halos that are significantly smaller than our estimate, we find that only two 
out of 16 objects with V max > 450 km sec -1 in the simulation have no sub-halos with 
M > 10 11 /i~ 1 M Q . We conclude that it would be very unlikely for an object like LBG-2377 
to exist without multiple components. 

By following the merger trees of our three V max > 600 km sec -1 halos, we find that 
they evolve from M = 1O 13 ' o_13 ' 5 /i _1 M0 halos at z = 3 to cluster-mass systems with M = 
1O 13 ' 9 ~ 14 ' 3 /i -1 M at z = 0. These are among the 12 most massive halos in the box at z = 
and contain ~ 10 — 50 sub-halos larger than lO n /i _1 M0. This analysis suggests that LBG- 
2377 will evolve into the dominant galaxy within a large group or cluster by the present 
day. 
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3.3. Ly-a emission 



The double peak nature of the Ly-a feature in LBG-2377 was only detected with higher 
resolution and SNR. Ho wever, multip l e Ly- a peaks are not unique to LBG-2377. In the LBG 
spectroscopic survey of ICooke et al.l (120051 ). we find ~ 3% of the low-resolution spectra to 
display double Ly-a emission peaks and resolved double peaks in their deep image profiles. 
We can only assume that this fraction increas es with higher resolut ion or SNR. One example 
of this effect may be found in the sample of IShapley et al.l (120061 ). Of the nine deep high 
SNR LBG spectra with Ly-a in emission, three have double peaks. Two of the remaining 
six have redshifts within Az = 0.01 and a separation of ~ 2" and are consistent with a 
double-peak system at a larger angular separation. Interestingly, upon inspection of the 
average internal velocity differences between the Ly-a emission and ISM absorption features 
(Az = z em — z a bs), the single peak systems are reported to have (Az) ~550 km sec -1 , whereas 
the double peak systems have (Az) ~ 1200 km sec -1 §. One plausible interpretation of this 
discrepancy is the result of the velocity separations, and consequent spectral blending, of two 
individual interacting systems. Confirming separate components in LBGs showing only Ly- 
a in absorption is more difficult. However, there are indications o f separate compone nts in 
the intermediate- resolution, high SNR spectrum of MS1512-cB58 (jPettini et al.ll2002l ). The 
prediction implied is that with higher SNR and/or resolution, a significant fraction, possibly 
as high as ~ 30%, of LBGs with Ly-a in emission have multiple massive components and 
may indeed be separate interacting or merging systems. This result would be consistent 
with the fraction of expected major m ergers at z ~ 3 from numerical simulations (e.g. 
Wechsler et al.l l200ll ; iRyan et al.l 120071 ) ). and with galaxy formation theories that suggest 



mergers are an important contributor to the star formation activity in LBGs (IKolatt et al. 
19991 : ISomerville. Primack fc Faberll200lh . 



4. Discussion 

What is the nature of LBG-2377? Possible interpretations are that LBG-2377 represents 
either (1) a single extended LBG with M ~ 10 13 M Q exhibiting multiple regions of vigorous 
star formation, (2) a single, extended LBG with mass more typical of LBGs (M ~ 1O 12 M ) 
that appears over-luminous due to exceptionally high bursts of star formation, (3) multiple 



3 This value omits the system C32 that shows a higher redshift for absorption than emission for one of 
the two components. Including this system results in a Az ~ 900 km sec -1 . If the absorption features of the 
lower redshift systems are weak or narrow in this interpretation, adopting a typical Az for the that system 
would result in Az ~ 1100 km sec -1 for the assumed two components. 
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unbound M ~ 1O 12 M LBGs in near proximity, or (4) multiple LBGs in a M ~ 1O 13 M 
halo observed in a late stage of merging. We discuss each scenario below and show that the 
evidence gathered to date favors interpretation 4. 

If the component stellar photospheric lines are misidentified and the complex nature of 
the ISM absorption and Ly-a emission lines are ascribed to violent mixing of the ISM and 
discrete star forming regions separated by > 2bh~ l kpc, interpretation 1 may hold. However, 
a sine; le M ~ 1O 13 M halo with no massive sub-halos is in conflict with simulation results. 
In this interpretation, LBG-2377 is indeed a massive and extended system that appears to be 
experiencing a monolithic burst of star formation. Under interpretation 2, LBG-2377 would 
be in less conflict with simulations yet the above assumptions for stellar and ISM lines must 
hold. In this case, LBG-2377 is a single system undergoing an unusually strong star formation 
episode at a level that has not been witnessed in any other LBG. In addition, this exceptional 
star formation must be occurring in multiple regions separated by > ThhT 1 kpc. If the 
velocity offsets are interpreted as redshifts, as in interpretation 3, the systems would have 
separations of ~ 5/i _1 Mpc comoving. An event of three or more unbound systems with M ~ 
1O 12 M , < l."0 angular separations, and an LBG density of ~ 10~ 3 — 10~ 4 /i 3 Mpc~3 has a 
vanishingly low probability, even if the systems were given reasonable peculiar velocities. The 
most likely scenario is interpretation 4 because it is consistent with the extended nature and 
multiple peaks seen in the images, the superposed spectra of multiple identified components, 
the observed velocity dispersion, and the expectations from the simulation. In addition, this 
interpretation provides a natural driver for the observed star formation and bright integrated 
luminosity. Furthermore, the double Ly-a emission peak of LBG-2377 is also seen in other 
LBGs which exhibit two distinct components in their deep images and Ly-a emission-line 
velocity differences consistent with close galaxy pairs. 

LBG-2377 is a fortuitous discovery that enables detailed study of the individual compo- 
nents of a candidate massive galaxy in the late stage of merging. Study of the relationships 
between the spectral and spatial properties of this single system will constrain the formation 
processes of the LBG population and lend insight into the fraction of LBGs detected as a 
result of merger induced star formation. Future high-resolution ground-based IFU data on 
existing instruments, such as OSIRIS on Keck, have the capability to resolve the velocity 
structure of LBG-2377 and hence its true nature. 
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for Cosmology at the University of California, Irvine. E.J.B. and J.S.B. are supported by 
the Center for Cosmology at the University of California, Irvine and J.S.B. and K.S. are 
supported by NSF grant AST-0507816. The authors wish to recognize and acknowledge the 
very significant cultural role and reverence that the summit of Mauna Kea has always had 
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Fig. 1. — Contour plots of LBG-2377. The 2-D plot indicates a primary peak located slightly 
to the E of the object centroid. Secondary and tertiary peaks are detected W and SE of 
the centroid, respectively. Three 3-D contour plots are inset with rotations (clockwise from 
upper-right), 185°, 335°, and 30° (S= 0°). 
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Fig. 2. — Spectrum of LBG-2377. This section is shown to illustrate both the double-peak 
Ly-a emission feature and broad, complex interstellar absorption features [vertical (blue) 
dashed lines]. In all sepctra in this paper, the green line indicates zero flux and the red line 
is the error array. 
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Fig. 3. — Spectrum of LBG-2377 showing five tentative components (A - E in Table [T]) 
identified using 21 ISM features. Velocities, are with respect to the average of components 
A, C, & D. Upper left: The two strong Ly-a emission lines (components C & D) and three 
potential weak Ly-a lines [A, B, & E]. Clockwise from the upper right: Fits for the ISM 
transitions Sill, Cn, and Ol. Also shown are the stellar photospheric features Sim 1294.5, 
1296.7, Cm 1296.3, Cn 1323.9, Nm 1324.3, and Oiv 1343.4 (not labeled for clarity) indicated 
by short vertical lines. The signal-to-noise ratio is ~ 15 in these regions, therefore the more 
prominant features are real. Each absorption feature is in agreement with their respective 
systemic redshift to within the spectral resolution (~ 150 km sec -1 ). 



-14- 



Table I. LBG-2377 Components 



Comp. 


ZPhot 


ZISM & 


ZLya 


A«/ 5 M a b 


Av Lya h 


A 


3.0385 


3.0354 


(3.0488) c 


-230 


(765) c 


B 


3.0343 


3.0308 




-260 




C 


3.0289 


3.0258 


3.0416 


-230 


945 


D 


3.0244 


3.0209 


3.0341 


-260 


725 


E 


(3.0161) c 


(3.0143) c 


(3.0263) c 


(-135) c 


(760) c 



a Averaged fit to multiple features. 

b Approximate velocity offset with respect to the stellar pho- 
tospheric velocities in km sec -1 . 

c Features are weak and respective values have a lower confi- 
dence level. 



